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Comment on "Non-Mean-Field Behavior of the 
Contact Process on Scale-Free Networks" 

Recently, Castellano and Pastor-Satorras [l[ used the 
finite size scaling (FSS) theory to analyze simulation 
data for the contact process (CP) on scale-free networks 
(SFNs) and claimed that its absorbing critical behavior is 
not consistent with the mean-field (MF) prediction. Fur- 
thermore, they pointed out large density fluctuations at 
highly connected vertices as a possible origin for non-MF 
critical behavior. In this Comment, we propose a scaling 
theory for relative density fluctuations in the spirit of the 
MF theory, which turns out to explain simulation data 
perfectly well. We also measure the value of the critical 
density decay exponent /8/fi|, which agrees well with the 
MF prediction. Our results strongly support that the CP 
on SFNs still exhibits a MF-type critical behavior. 

In contrast to equilibrium models with thermal noise, 
nonequilibrium models exhibiting absorbing critical phe- 
nomena involve so-called multiplicative noise. In partic- 
ular, the noise amplitude is proportional to the square 
root of the particle density in models belonging to the 
directed percolation class such as CP [2]. Therefore, at 
the MF level neglecting spatial fluctuations, the particle 
number fluctuation should be proportional to the square 
root of the total number of particles. It is easy to derive 
the relative density fluctuation on vertices of degree k as 

r k = A Pk / Pk ~ [NpkPk]- 1 / 2 , (1) 

where p k is the partial density on vertices of degree k, N 
is the total number of vertices, and Pk is the degree dis- 
tribution. As pk ~ pk near criticality (p: total density), 
we find, for reasonably large k, 

r k ~ {Np)-^k^'\ (2) 

where 7 is the degree exponent. 

For large enough k > N 1 ^ 1 where the degree distribu- 
tion becomes almost flat in practice (NPk ~ 0(1)), we 



FIG. 1: (Color online) Data collapse of the relative density 
fluctuations rk averaged over surviving samples in the steady 
state for 7 = 2.75 (p c = 0.4240, c = 0.363, y = 0.104) and 
7 = 4.0 (p c = 0.3570, c = 0.25, y = 0.125), respectively. 
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FIG. 2: (Color online) Effective exponents f)/v\\ for 7 = 2.75 
with N = 10 7 . All samples are alive when data are obtained. 

expect rk ~ (p/c) -1 / 2 . With these two limiting results 
and also with p ~ N~ a at criticality, one may assume a 
critical scaling relation as 

r k =iV (Q - 1/7)/2 /(fc^ _1/7 ), (3) 

where f(x) ~ xlT -1 )/ 2 as x — > and f(x) ~ a; -1 / 2 as 
x — > 00. In Fig. [1] our simulation results at criticality 
for 7 = 2.75 and 4.0 with various sizes N — 10 3 , • • • , 10 6 
show excellent agreements with our prediction. Note that 
the case of 7 = 4.0 where the 7-independent ordinary MF 
appears also shows the similar behavior in rk- 

We also observe the density decay dynamics for 7 = 
2.75 in Fig. [2j where the effective decay exponents /3/v\\ 
are plotted against 1/t. Our best estimate is /3/Vii = 
1.34(5) which agrees well with the MF prediction of v\\ = 
1 and /3 = 1/(7-2) for 7 < 3. 

Finally, we comment on the FSS exponent a = f3 jv. 
Vertices of degree k > iV 1 / 7 are practically irrelevant. 
Accordingly, the cutoff in k bigger than iV 1 / 7 can not 
influence critical FSS which should be determined by 
the competition between correlated volume and system 
size. We conjecture that v = (7 — l)/(7 — 2) from the 
hyperscaling-type argument Q. Our numerical estimate 
j3Jv — 0.59(2) for 7 = 2.75 agrees well with our conjec- 
ture. In fact, the results of Ref. [lj also agree reasonably 
well with our conjecture. In conclusion, we believe that 
the CP on SFNs exhibits the MF-type critical behavior. 
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